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Abstract 



' Motivated by a recent experiment by Weiss et al. [Phys. Rev. Lett. 70, 4118 



' (1993)], we present a detailed study of quantum transport in large antidot 

arrays whose classical dynamics is chaotic. We calculate the longitudinal and 
Hall conductivities semiclassically starting from the Kubo formula. The lead- 
ing contribution reproduces the classical conductivity. In addition, we find 
oscillatory quantum corrections to the classical conductivity which are given 
' in terms of the periodic orbits of the system. These periodic-orbit contri- 

butions provide a consistent explanation of the quantum oscillations in the 
magnetoconductivity observed by Weiss et al. We find that the phase of the 
oscillations with Fermi energy and magnetic field is given by the classical ac- 



^ ' tion of the periodic orbit. The amplitude is determined by the stability and 



the velocity correlations of the orbit. The amplitude also decreases exponen- 
tially with temperature on the scale of the inverse orbit traversal time fi/T^. 
The Zeeman splitting leads to beating of the amplitude with magnetic field. 
We also present an analogous semiclassical derivation of Shubnikov-de Haas 
oscillations where the corresponding classical motion is integrable. We show 
that the quantum oscillations in antidot lattices and the Shubnikov-de Haas 
oscillations are closely related. Observation of both effects requires that the 
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elastic and inelastic scattering lengths be larger than the lengths of the rel- 
evant periodic orbits. The amplitude of the quantum oscillations in antidot 
lattices is of a higher power in Planck's constant h and hence smaller than 
that of Shubnikov-de Haas oscillations. In this sense, the quantum oscillations 

in the conductivity are a sensitive probe of chaos. 
73.50.Jt, 03.65.Sq, 05.45.+b 
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Typeset using REVT^ 



I. INTRODUCTION 



Over the past few years high-mobility semiconductor heterostructures have proved to be 
a valuable tool for studying the implications of quantum chaos experimentally [|l|. Advances 
in technology have made it possible to manufacture structures whose elastic and inelastic 
mean free paths exceed the relevant sample dimensions. In these ballistic (micro) structures 
essentially all large-angle scattering is due to the geometric features of the device . Hence, 
they can be used in particular for experimental studies of billiard systems which have 
been central to the theoretical development of quantum chaos [|I| . These novel experimental 
realizations of quantum- chaotic systems are important because there are actually only a 
small number of systems in the field of quantum chaos that allow for a controlled comparison 
between theory and experiment. 

Much of the work on quantum chaos has focused on the spectral properties of systems 
whose classical dynamics is chaotic 0]. However, the detailed spectrum is not an accessi- 
ble observable in ballistic semiconductor structures. Rather, one measures thermodynamic 
properties such as the magnetic susceptibility and the persistent current or transport 
properties such as the conductance of the sample. Hence the challenge to theory was 
whether these quantities are sensitive to the nature of the underlying classical dynamics. 
This question has now been answered in the affirmative for both thermodynamic and trans- 
port properties. With regard to transport properties of ballistic microstructures, it has been 



predicted that the power spectrum of the conductance fluctuations |110| , p!l| and the lineshape 



of the weak-localization peak [TllJT^ are probes of quantum chaos and there is accumulat- 
ing experimental evidence in favor of these predictions 0,^,0. There are also predictions 
that the amplitudes of the persistent current |jl3| and the magnetization |ll4| , p^ are sensitive 
probes of quantum chaos. These have not yet been subject to experimental tests. 

Chaotic transport in semiconductor structures has been investigated in two kinds of 
experiments. Experiments on mesoscopic ballistic cavities with the shape of chaotic billiards 
coupled to leads HflM] have focused on generic effects such as conductance fluctuations 



3 



and weak localization. These effects are familiar from the study of transport in diffusive 
conductors [|T6|-|T^. In the other type of experiment one measures magnetotransport in 



macroscopic antidot lattices. Antidot lattices consist of a periodic array of holes "drilled" 
into a two-dimensional electron gas P,|5|JT9||. They are macroscopic in the sense that the 



phase-coherence length is much smaller than the system size (the lattice spacing of the 
antidot array, on the other hand, may be smaller than the phase- coherence length) 
Antidot lattices can be viewed as an approximate experimental realization of the Sinai 
billiard |^0[| which is a well-known chaotic system. Experiments on antidot arrays were 



mostly concerned with features in the magnetotransport which reflect the specific nature of 



the sample. In this paper we present a theoretical study of transport in antidot lattices p2 

Two analytical methods dominated theoretical studies of quantum chaos, random-matrix 
theory and the semiclassical approximation . While the first is most appropriate for 



describing generic effects due to long trajectories, semiclassical theory is especially useful in 
describing effects due to short system-specific trajectories. We find that semiclassical theory 
is particularly suited to address the system-specific features in the magnetoconductivity of 
antidot lattices. In the context of quantum chaos interest in semiclassical theory derives 
mostly from the fact that, strictly speaking, chaos is a concept of classical physics. Semi- 
classical theory allows one to study most explicitly the implications of (classical) chaos in 
quantum mechanics. In the context of ballistic structures semiclassical theory is a natural 
approach because the Fermi wavelength is often the shortest length scale in the problem. 
Furthermore, the contributions of long classical trajectories are frequently suppressed by 
finite temperature or inelastic scattering, rendering the semiclassical expressions quite man- 
ageable [|13|-[T5|1. 



The first experiments on antidot lattices by Weiss et al. focused on the regime of 
classical transport. They revealed a series of pronounced magnetoresistance peaks at mag- 
netic fields where the classical cyclotron radius is commensurate with the lattice spacing of 
the antidot array. For these magnetic-field values, the electrons can circle around a certain 
number of antidots for long times. This effectively slows them down leading to an increase 



in the resistivity P J25|J26[] . Here we are concerned with more recent experiments at lower 
temperatures |^ which revealed additional quantum oscillations in the magnetoresistance 
superimposed on the classical peaks. The period of these oscillations corresponds to roughly 
a flux quantum threading a unit cell of the antidot lattice. Depending on the sample these 
oscillations varied between 1/ B periodicity familiar from Shubnikov-de Haas oscillations 



and B periodicity familiar from Aharonov-Bohm oscillations Weiss et al. explained 



these observations phenomenologically by assuming that the resistivity is proportional to 
the square of the density of states, pxx ~ d^iEp, B). The period of the oscillations can then 
be understood in terms of a few short periodic orbits which, according to Gutzwiller's trace 
formula p4|, give oscillatory contributions to the density of states. This approach is clearly 



unsatisfactory because it does not make any predictions concerning the absolute amplitude 
of the oscillations and, more importantly, because this relation between the resistivity and 
the density of states cannot be derived from first principles. In fact, this relation certainly 
does not reproduce the pronounced classical peaks in the resistivity p9| . 

Nevertheless, the apparent success of the phenomenological theory of Ref. ^ suggested 
that there is indeed a relation between quantum corrections to transport properties and 
classical periodic orbits. In this paper we present a semiclassical theory for the conductiv- 



ity starting from first principles |^. Since we are considering macroscopic antidot lattices 
much larger than elastic and inelastic mean free paths we start from the Kubo formula. To 
leading order in Planck's constant h we reproduce the classical Kubo formula. According 
to Fleischmann et al. p5| this contribution accounts for the classical peaks in the resistiv- 
ity. The central result of this paper is that we find quantum corrections to the classical 
conductivity which explain the additional quantum oscillations observed by Weiss et al. 
0. These quantum contributions to the conductivity are indeed given in terms of classical 
periodic orbits. The contribution of each orbit oscillates as function of Fermi energy and 
magnetic field with a phase equal to its classical action. The amplitude depends on the 
stability of and velocity correlations along the orbit. The predictions of our theory for the 
period, the amplitude, and the temperature dependence of the periodic-orbit contributions 



are consistent with the experimental results of Weiss et al. 0. It was previously known that 
periodic orbits are important in describing the oscillatory corrections to the density of states 
and hence for thermodynamic properties of mesoscopic systems [ pl| , p!^ - p^ , |32| , |33[] . Our 



results combined with the experimental results of Ref. show that periodic orbits are also 
important for quantum corrections of transport properties MM. There have also been at- 



tempts to interpret the "frequency spectrum" of conductance fluctuations in billiard-shaped 
ballistic microstructures in terms of periodic orbits 0,^. The relevance of our results to 
these experiments is not clear because the conductance of such phase-coherent samples is 
more appropriately described using the Landauer-Biittiker approach rather than the Kubo 
formula. 

We find that the quantum oscillations in the conductivity of "chaotic" antidot lattices are 
closely related to the Shubnikov-de Haas oscillations in the conductivity of the "classically 
integrable" two-dimensional electron gas. To make this relation explicit we present an 
analogous semiclassical derivation of Shubnikov-de Haas oscillations. Even though both 
effects are closely related, the technical treatment is somewhat different. It is well known 
that integrable and chaotic systems must be treated differently in the semiclassical theory 
of the spectral density Likewise we find that certain spatial integrals which must 



be computed exactly in the derivation of Shubnikov-de Haas oscillations are evaluated by 
the stationary-phase method for chaotic antidot lattices. As a result, the amplitude of 
Shubnikov-de Haas oscillations is larger than that of the quantum oscillatons in antidot 
lattices. In analogy to the semiclassical theory of the spectral density [|2^, our results for 



antidot lattices are strictly valid only if the classical dynamics is completely chaotic. In 
section [III C] we present evidence that this is indeed the case to a good approximation for 
the main sample investigated in Ref. [Q. 

The outline of the paper is as follows. In section |1| we present the semiclassical theory for 
Shubnikov-de Haas oscillations. Section |T| contains the central results of this paper on the 
longitudinal conductivity ( [111 A| ) and the Hall conductivity ( [111 B| ) of antidot lattices, and 



a comparison between our theory and the experimental results of Ref. |^ ( [HI Q) . Various 



details of the calculations in sections |T|, [III A| , and [III B| are given in three appendices. 
Finally, we conclude in section 



II. SHUBNIKOV-DE HAAS OSCILLATIONS 

In this section we present a semiclassical derivation of the Shubnikov-de Haas oscillations 
in the magnetoconductivity of two-dimensional electron gases. Our motivation for doing so 
is twofold. First, it is instructive to see how the semiclassical methods which we employ 
in our derivation of the quantum oscillations in antidot lattices work in a more familiar 
context. Second, it turns out that the Shubnikov-de Haas and the quantum oscillations in 
antidot lattices are closely related phenomena. This relationship is brought out explicitly 
by comparing the calculations for both effects. 

Throughout this paper we consider transport in macroscopic systems in the sense that 
the sample dimensions are much larger than the elastic and inelastic mean free paths. In 
this regime, transport is appropriately described by the Kubo formula. Here, we start from 
the Kubo formula for the longitudinal conductivity 



■ Tr {vJriEp - H)vMEp - H)} , (1) 



n 

where H = (l/2m)(p — eA(f))^ is the Hamiltonian of a two-dimensional electron gas in 
a perpendicular magnetic field, A = {—By/2, Bx/2) is the vector potential taken in the 
symmetric gauge, and is the x component of the velocity operator (here m denotes the 
effective mass of the electrons). The Fermi energy is denoted by Ep and the volume of the 
system by f2. It is well known that one needs to include some disorder scattering when 



directly computing the dc conductivity within the Kubo formalism |3^. We include weak 
disorder at the level of the Born approximation by giving the 5-functions a finite width 
r = h/2Tei {Tel is the elastic scattering time). This is sufficient to prove our main point, i.e., 
that there is a precise connection between classical periodic orbits and quantum corrections 
to the conductivity. However, it should be kept in mind that the Born approximation cannot 
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fully describe smooth disorder potentials because it does not distinguish between the elastic 
and the transport mean free paths. Unfortunately, it appears to be a difficult problem to 
go beyond the Born approximation when using the semiclassical approach presented in this 
paper. In principle, the level broadening F also includes a contribution due to inelastic 
scattering. This will be left implicit throughout the paper. 

To evaluate the Kubo formula we start by expressing the 5-functions in Eq. (D) in 
terms of the retarded and advanced Green functions G^{E) = [E ~ H ± *r]~^, 

1 



SriEp - H) 



2ni 



[G+{Ef)-G-{Ef)]. 



(2) 



Evaluating the trace in ([^) in position representation, one finds that the longitudinal con- 
ductivity can be expressed as 



^xx ^xx ~^ ^xx ) 



(3) 



where 



e^nh/ 1 ^2 



\2nmJ 



-V, 



\Gf ,{E]: 



and 



e^Tch/ 1 



Q \27im 



'-V^-eA^{v)\G-,,{EF) 



-V^-eA,{v')\Gt,(EF) 



(4) 



V'^-eA,{v')^Gt,^AEF) 



(5) 



Here, Gfj.,{E) = {r\G^(E)\r') denotes the matrix elements of the Green functions. It is 
shown in appendix y that both a^^ and vanish in the semiclassical approximation. 
Hence, we focus on a^^ in the remainder of this paper. 

Having expressed the longitudinal conductivity in terms of Green functions, we now 
approximate these by their semiclassical expressions. Generally, the semiclassical Green 
function G^, ^{E) is a sum over all classical trajectories of energy E from r to r'. For free 
electrons in a magnetic field, there are an infinite number of trajectories connecting r to r' 
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if the distance between r and r' is less than twice the cyclotron radius Re- Two different 
cyclotron orbits labeled by S (for short) and L (for long) which can each be traversed 
completely for n = 0,1,2... times pass through both points (cf. Fig. 1). There are no 
trajectories from r to r' when |r — r'| > 2Rc and the semiclassical Green function vanishes 
for such arguments. An explicit calculation shown in appendix ^ gives 

Here Tn,q denotes the traversal time of the trajectory (n, g). 



21111 , . 

^n,q I" tqi ( ' ) 



where 



2 

ts = — arcsin(|r - r'\/2R^) 
2 

= — [tt - arcsin(|r - r'|/2i?c)] (8) 

are the times for going from r to r' without completely traversing the cyclotron orbit (10^ = 
eB/m is the cyclotron frequency). Each trajectory contributes with a phase factor given by 
its classical action 



Sn,q{E) = ETn^q H r""^ 



r — r f UrT, 



cot ( ^""^-1 ^ _ {^xy' - x'y) 



(9) 



and the Maslov indices rjn^s = 2n, ^ = 2n + 1. The corresponding expression for G^. j.i{E) 
follows from the relation 

G-AE) = [Gt,AE)r. (10) 

One notes that both Green functions appearing in (^ are given in terms of trajectories from 
r to r'. 

Substituting these semiclassical expressions for the Green functions into Eq. (|^) for the 
conductivity one has a double sum over classical trajectories [n, q) from r to r'. 
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2nE 



ni=0 n2=0 qi=S,L q2=S,L 



where we have used the relations 



m \ I 



I (n 



m \i 



(12) 
(13) 



Here {vq)x and {v'^x denote the x components of the velocities at the points r and r', 
respectively. One notes that for qi = q2, the difference in action Sn2,q2 —Sni,qi is independent 
of both r and r'. Only these terms contribute to leading order in the semiclassical limit 
because the integrand oscillates rapidly when q2- A further simplification arises because 
both the times T„,g and the difference in action depend only on |r— r'|. It is useful to introduce 
center-of-mass coordinates and to change variables to the times tq, 



[ (fr f d\' = 27rn / ^ |r - r'l d\r - r'\ 

Jn J\r-r'\<2Rc Jo 



/ dts\sm{LOcts)\+ i (iti I sin (c^c^l) 



(14) 



Using {yq)x{y'g)x = Vx{S^)vx{tq) and combining the integrations over ts and t^, one obtains 



27ih^ Jo 

00 

X 

ni, 712=0 



dtvx{0)vx{t) e-*/2^^' 



00 f 



-(ni + n2)| cos |27r - (ns - ni)| . (15) 



Finally performing the sum over rii + and introducing p = \ni — we find 



2e2 / Ep 



h V 2'KhLJr 



C{Vx,Vx 



1 + 2 y cxp 



cos < Znp - 



(16) 



Here we introduced the correlation function 



C{Vx,Vx) 



1 f2iT/oJc roo 

/ dt / dt' Vx(t)vx(t + t') e-^ . 
Jo Jo 



RlJo 



(17) 
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This way of writing the result will be useful for comparison with the conductivity of the 
antidot lattice derived in the next section. For comparison with the standard result for 
Shubnikov-de Haas oscillations in the literature |^ one easily computes 

C{v^,v^) = —- -. 18) 

1 + [uJcTelY 

Using that the density of electrons (per spin) is Ue = mEp/2iTh'^ the prefactor in Eq. (p!6D 
can be rewritten in the more familiar form (?2ee^re;/m)[l + (ucTei)'^]^'^ . This proves that 
the semiclassical approximation reproduces the low-field approximation {Ep ^ huc) to the 



quantum-mechanical result (within the Born approximation) |2^. We note that the classical 
Drude conductivity arises from the diagonal terms in the double sum over trajectories in 
Eq. (pJ]). The Shubnikov-de Haas oscillations arise from non-diagonal terms corresponding 
to different numbers rii and n2 of complete revolutions around the same cyclotron orbit. In 
the next section we will find a closely analogous structure of the terms contributing to the 
conductivity of antidot lattices. 



III. QUANTUM TRANSPORT IN ANTIDOT LATTICES 

A. Longitudinal conductivity 

In the following we consider the longitudinal conductivity of macroscopic antidot lattices. 
We start from the Kubo formula (|l]) with the Hamiltonian H = (l/2m)(p — eA)^ + V^(r) 
where V"(r) is the antidot potential. A semiclassical approximation to the Green function 
G^, j.{E) for general potentials V{r) is obtained by making consistent use of stationary-phase 
integration both for the path integral for the propagator and for the Fourier transform to 
the energy domain. One finds that the Green-function matrix elements are approximated 
by sums over the classical trajectories 7(r, r') of energy E going from r to r' ||24|| , 

GJ.W = 7^(^)'" E fl,(-.r>xp{is,(r.r';£)-»^„,}. (19) 
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The phase of the Green function is determined by the classical action Sj{r,r'; E) = J^drp 
of the trajectory 7 (p is the canonical momentum and the integral is taken along the tra- 
jectory) and by the Maslov index rj^, which counts the number of conjugate points along 
the trajectory. Each trajectory contributes with an amplitude which is determined by 
disorder and the classical phase-space density, 



1/2 



D^{r,r') = exp y-^^^ 



dr'-drj dEdrj 



(20) 



Here the traversal time of the trajectory is denoted by T^. The corresponding expression for 
Gj^j.,{E) follows from G~j.,{E) = [G^, ^^.{E)]* . Importantly, both Green functions appearing in 
(§) are given in terms of trajectories running from r to r'. The main difference between the 
semiclassical expressions for G^, j.{E) and G^^y^E) is that the actions enter in the exponent 
with opposite signs. 

Inserting the semiclassical expressions for the Green functions into Eq. (|^) one finds 
that (Jxx is given by a double sum over classical trajectories from r to r'. In the simplest 
approximation one retains only the diagonal terms in this double sum for which the rapidly- 
oscillating phase factors drop out, 

,2_fc / 1 X 3 



' XX — 



^^(i) I 'PrSr' j: .^KD^;(rX). (21) 

7(r,r') 



Here we used the elementary relations 

- ( -V. - eA,{r)] S^{r, r'; Ep) = -v, (22) 

^ (^V; - eA^ir')^ 5,(r, r'; Ej.) = v'^, (23) 

where and v'^ denote the x components of the initial and final velocities, respectively. 
This contribution gives the classical conductivity. To see this, we multiply the right-hand 
side of Eq. ( |2TD by 1 = / dE 6{E — Ep) and introduce as new integration variables the 
initial position r, the initial momentum p, and the traversal time t of the trajectory 7. The 
corresponding Jacobian is given by 
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d?rd^pdt = d?rd^r' dE ^ 

7(r,r') 



dpi dpi 
dt dt 



(24) 



The sum over 7 appears because initial conditions (r,p) specify a classical trajectory uniquely 
while boundary conditions (r, r') can be satisfied by more than one trajectory. One notes 
that the Jacobian (p^) involves precisely the same determinant which also appears in the 
amplitude factors D^{r, r') in the Green functions ([T9|). Hence, one finds for the conductivity 

-,2 



XX 



^ j d\ d^pdtvMvxit) e-"^^' 5{Ef - H{^,v)) , (25) 

where i?(p, r) denotes the classical Hamiltonian. Introducing the classical phase-space av- 
erages 

(/(P,r))r = j;^^^ J d'rd'pf{p,r)S{E^-H{p,r)), (26) 

(A^(0) = m/2Tih^ is the density of states at the Fermi energy in two dimensions) the con- 
ductivity can be rewritten in the familiar form 



POO 

e^N{0) dt{vMvx{t))ve-"^^K (27) 



This is precisely the classical Kubo formula. Hence, the diagonal approximation to the 
double sum over classical trajectories in Eq. (|^) is equivalent to the classical approximation. 
At first sight, it may appear inconsistent that the classical limit (|27| ) still - via the density of 
states A^(0) - involves Planck's constant h. However, this is simply a consequence of the fact 
that we use the Fermi-Dirac distribution instead of the Boltzmann distribution appropriate 
for a complete classical limit. 

As shown by Fleischmann et al. p5| the classical Kubo formula ( pTl) explains the peaks 
in the magnetoresistance when the cyclotron radius and the lattice spacing of the antidots 
are commensurate. In this paper, we are mainly interested in the quantum corrections 5(Jxx 
to the classical conductivity. Clearly, they must arise from the nondiagonal terms in the 
double sum in (^, 
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7l(r,r')7^72(r,r') 

X exp 1^ [S^^{r, r'; Ep) - S^,{r, r'; Ep)] - y [r]^^ - ^71]} • (28) 



In these terms, the phase factors do not cancel and hence the integrand involves rapidly- 
oscillating factors. Within the semiclassical approach for completely chaotic systems, it is 
natural to evaluate the spatial integrals by the stationary-phase method. This procedure 
becomes exact in the semiclassical limit h ^ 0. Considering the contribution arising from 
trajectories 71 (r, r') and 72 (r, r'), the stationary-phase conditions for the r and r' integrations 
are 

V[S,,{r, r'; Ep) - S,,{r, r'; Ep)] = p, - = (29) 
V'[5^,(r, r'; Ep) - S,,{r, r'; Ep)] = p', - p[ = 0. (30) 

Here, Pi and p- denote the initial and final (canonical) momenta of trajectory 7j. Since the 
vector potential cancels from the stationary-phase conditions, one finds that the integrand 
becomes stationary when initial and final velocities of 71 and 72 are identical. Since also the 
initial and final positions of the two trajectories must each coincide one might at first think 
that this implies 71 = 72 which would just be the classical contribution discussed above. 
However, there is an additional possibility to satisfy the stationary-phase conditions ( PP||5D[ ) 



when the two trajectories are part of a classical periodic orbit. Then 71 and 72 can differ 
by an integer number p of traversals of the orbit as illustrated in Fig. 2. In fact, for each 
(primitive) periodic orbit 7 there is an infinite set of pairs of trajectories 71 and 72 which 



satisfy the stationary-phase condition (|29| , [30| ). The pairs are labeled by p which was defined 
above and by n which counts the number of times the shorter trajectory from r to r' returns 
to its starting point r. It is these quantum corrections to the conductivity which we focus 
on in the following. 

The evaluation of the stationary-phase integrations over r and r', while different in 
detail, is closely analogous to performing the final trace in the derivation of Gutzwiller's 
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trace formula ||2^. First we note that the saddle-point action is just given by the action 
S^{Ep) of the (primitive) periodic orbit, 

S^,{r, r'; Ep) - S^,{r, r'; Ep) = pS^iEp). (31) 

To evaluate the integral over the quadratic fluctuations around the stationary-phase point 
it is useful to introduce a new set of orthogonal coordinates r = {z, y) where z measures the 
arclength along the periodic orbit 7 and y measures the transverse deviations from the orbit 
(cf. Ref. [^]). The Jacobian of this transformation is unity, d'^rd'^r' = dzdydz' dy'. While 
the integrals over y and y' can be done by the stationary-phase method, the integrals along 
the periodic orbit must be performed exactly because the second variation of the action 
vanishes in this direction. Introducing the matrix F of second derivatives of S'21 = 5*^2 ~ "^-n 
with respect to y and y', 



(32) 



d^S2l dydy,S 21 

dyidyS2i dy,S2i 

(here we used the shorthand notation dy = d/dy) and the vector = [y,y'] of transverse 
deviations from the periodic orbit, one finds for the fluctuation integral 

exp{A,.fy}^-^. (33) 

A straight-forward but tedious calculation presented in appendix |B] shows that this fluctu- 
ation amplitude combined with the semiclassical amplitudes of the Green functions can be 
expressed in terms of the monodromy matrix of the periodic orbit alone, 



(34) 



detF| |v(2)||v(z') 

Finally, we simplify the integrations along the periodic orbit over z and z' . Since the 
integrations over r and r' are over the volume of the system, the integrations over z and z' 
extend over a full period around the primitive periodic orbit. Furthermore, since r is the 
starting and r' the end point of the trajectories, the limits of integration are < z < 
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and z < z' < z' + L^, where is the length of the (primitive) periodic orbit. Hence, the 
z' integrations become 

'+^^ , / v.^{z)v,^{z') 



= exp{-pT^/2Tei) V / dt dt'v^{t)v^{t + t') exp{-{n^ + t')/^e«} 



n=0 



pl'y poo 

= exp{-pT^/2T,i) dt dt'vS)vS + t')exp{-t'/Tei). (35) 
JO Jo 

Here, = L^/vp denotes the period of the primitive periodic orbit. Due to the translational 

symmetry of the antidot lattice each periodic orbit appears Q/a'^ times (a is the lattice 

spacing). Including only a single representative of each class of periodic orbits in the sum 

over 7, we find the final result 

2e2i?2 os^ cos{pS^{EF)/n - 'Kpa^/2} 
Sa^x = -f^-^2^2^ exp{-pT^/2Tei)C^{v^, f J |TrM^ - 2|V2 ' ^^^^ 

where the sum over 7 is over all primitive periodic orbits and p counts the repeated traversals. 
We also introduced the correlation function 

C^{vx,v,) = — dt dt'v,{t)v^{t + t')exp{-t'/Tei). (37) 
Jo Jo 

Some periodic orbits of antidot lattices turn out to be quite close to circular. For these it 
can be a useful approximation to replace the present correlation function by its analog {\[7\] 
for Shubnikov-de Haas oscillations. By collecting all phases one convinces oneself that the 
Maslov index appearing in Eq. (|36|) is the same that enters Gutzwiller's trace formula 



2J]. Hence, periodic-orbit contributions to the conductivity and to the density of states 
have the same phase which is physically reasonable. The result (|36D for the periodic-orbit 
corrections to the conductivity of antidot lattices should be compared to the result (p!6D 
for Shubnikov-de Haas oscillations. The principal difference is that the stability of the 
periodic orbit expressed in terms of the monodromy matrix enters in the case of chaotic 
dynamics. Furthermore, one notes that the quantum oscillations for the antidot lattices are 
of a higher power in Planck's constant h and hence smaller in the semiclassical limit than 
the Shubnikov-de Haas oscillations. 
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At finite temperatures the Kubo formula takes on the form 



n 



r ^^dE^^ {^^'^^^^ " H)vJr{E -H)}, (38) 



where /i denotes the chemical potential and ffj,{E) is the Fermi-Dirac distribution f^{E) — 
[1 + exp{P{E — fJ,))]~^ {(3 = l/ksT is the inverse temperature). To compute the temper- 
ature dependence of the periodic-orbit contribution to the conductivity one notes that the 
dominant energy dependence (to leading order in h) is in the phase factor. Hence, we need 
to compute the integral 

I{T) = - f dE cos [^S,{E) - ^pa,] . (39) 

Since the derivative of the Fermi function is strongly peaked at the Fermi energy, we extend 
the lower limit of integration to — oo and expand the action to first order around the Fermi 
energy, 

I{T) = Re— exp |-p5^(/x) - -pa^^ dE f^{E) exp [jVT^{E - /x)| . (40) 

Here we used that dS^/dE — T^. This integral can be evaluated by contour integration. 
The Fermi-Dirac distribution has poles at E — + i{2n + l)7r//3 with residues — One 
finds 

I{T) = cos[-SM - -pa,\ ^^^^^^^^i^^y (41) 
Thus we find for the periodic-orbit contribution to the conductivity at finite temperature 

^c^xx = -r — iJ2Yl exp{-pT^/2Tei)C^{v:c, Vx) 

{npT^/hp) cos{pS^{EF)/h - 7rpaj/2} 
smhiirpT^np) |TrM^-2|V2 ■ ^ ) 

The contribution from each periodic orbit decreases exponentially with temperature on the 
scale of h/Tj. In principle, there is an addional temperature dependence due to the variation 
of the inelastic scattering length with temperature. However, this temperature dependence 
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is slow compared to that arising from thermal smearing. For this reason, Eq. p^) does 
indeed describe the dominant temperature dependence of the effect. 

The Zeeman coupling of the electronic spin to the magnetic field is readily incorporated 
into our results. It effectively leads to separate Fermi energies Ep ± yUs-B for the two spin 
components (here /i_B = e%/2mf. is the Bohr magneton with mg the bare electron mass). 
Expanding the action to linear order in the Zeeman energy one finds that it leads to beating 
of the oscillation amplitude, 

(^CTxx = -f^^^Y^ exp(-pT^/2re;)C^(t;^, v^) 

The corresponding magnetic- field scale is Ti/ usT-y. 

B. Hall conductivity 

The calculation for the Hall conductivity a^y is closely analogous to that for the longitu- 
dinal conductivity described in the previous section. The Hall conductivity can be expressed 



in terms of two contributions 36 



(^xy = + a^J, (44) 



where 



< = ^Tr {v.,G^{EF)Vyb^{Ep - H) - vJriEp - H)vyG-{EF)} (45) 



2n 



and 



II 

^xy 



Here, p{E) denotes the density of states. First, we evaluate the contribution a^y using 
Gutzwiller's trace formula which gives a semiclassical approximation to the density of states 
p{E). The Weyl contribution to the density of states is independent of the magnetic field. 



Thus, we only need the oscillatory contributions Sp{E) to the density of states which can 
be expressed in terms of a sum over the periodic orbits 7 of the system, 

= ^^^g |TrM,^-2|V2 cos|-^^,(E) - . (47) 

Noting that to leading order in Planck's constant h we only need to keep the magnetic-field 
and energy dependence in the oscillatory factor, one readily finds 

We now proceed to evaluate the first contribution a^y. The evaluation of this term is 
completely analogous to that for . We do not consider the classical Hall conductivity 
and focus on the periodic-orbit contribution right away. Again, we only need to consider 
contributions from products of advanced and retarded Green functions. By analogy with 
the evaluation of a^y' we can immediately write down the result of the calculation for the 
Hall conductivity, 

^^xy = L exp(-pT^/2reO C^{v,, Vy) |TrM^-2|V2 ' (^9) 

where the velocity correlations are given by 

C^(^^x,t^j;) = -^y^ dtj^ dt'v^{t)vy{t + t')e~''^^K (50) 

The Maslov index is again the same that also appears in Gutzwiller's trace formula. 
Finite temperature and the Zeeman splitting lead to the same additional factors which also 
appeared in the treatment of the longitudinal conductivity in the preceding section. For an 
approximately circular periodic orbit of radius R traversed at constant angular velocity u) 
one finds 

^""^y^ h 1 + (cuTeO^ """P^ ^^''^^^'^ |TrM^-2|V2 " ^^^> 

The same disorder-dependent prefactor appears in the theory of Shubnikov-de Haas oscilla- 
tions [^]. Note that the amplitude of a circular orbit in the Hall conductivity is smaller by 
a factor 1/ujTf.i than that in the longitudinal conductivity. 
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C. Comparison with experiment 



Weiss et al. PJ37[] measured both the longitudinal conductivity and the Hall conductivity 
of macroscopic antidot lattices as function of magnetic field. The antidot lattices were 
patterned from high-mobility GaAs-AlGaAs heterostructures with elastic and inelastic mean 
free paths of order 5/im. For the main sample discussed below the lattice spacing was 
a = 200nm and the antidot diameter d = lOOnm. The Fermi wavelength was of order 
Xp = 50nm. 

At higher temperatures (typically T = 4.7K) Weiss et al. |^ found a number of peaks in 
the longitudinal resistivity at magnetic fields where the classical cyclotron radius is commen- 
surate with the lattice spacing of the antidot array. This phenomenon has been discussed 



theoretically by Fleischmann et al. ||2^. Here we are concerned with the additional quan- 
tum oscillations in both longitudinal and Hall conductivities superimposed on the classical 
conductivity. These were observed by Weiss et al. at subkelvin temperatures. The ex- 
perimental results for the main sample investigated in Ref. are shown in Fig. 3. To focus 
on the quantum oscillations, we have plotted 6a = a{T = 0.4K) — a(T = 4.7K). In the 
following we compare our theoretical results to the experimental data for this particular 
sample. The quantum oscillations have a period of the order of one flux quantum per unit 
cell and the oscillations in a^x and a^y are in phase. One reads off from Fig. 3 that the 
amplitude is roughly 6axx ~ 6axy ~ 0.015(kfi)~^. 

A detailed comparison between theory and experiment is complicated by a number of 
factors. Corrections to the semiclassical results may be relevant because the Fermi wave- 
length is not much smaller than the antidot radius. In addition, our approximations are 
strictly valid only for completely chaotic systems. By contrast, the experimental system 
exhibits mixed dynamics in part of the magnetic-field range over which the quantum oscilla- 
tions have been observed. One expects that this does not affect the phase of the oscillations 
strongly but does have an influence on their amplitude. A detailed comparison of theory 
and experiment also requires use of a model for the antidot potential which is not known a 
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priori. Perhaps the most serious shortcoming is our treatment of disorder which is sufficient 
to prove that there is a rigorous connection between quantum corrections to the conductivity 
and classical periodic orbits but which is not fully adequate to describe the smooth disorder 
potentials typical of semiconductor heterostructures. For want of a better treatment we have 
considered the elastic mean free path as a fit parameter, keeping it merely within reasonable 
bounds. 

We model the antidot potential by ^5 



V{r) = Vo [cos(7rx/a) cos{TTy/a)f . (52) 

Here, Vq controls the diameter of the antidots at the Fermi energy and /3 fixes the steepness 
of the antidot potential. Hence, all parameters but (3 are fixed directly by experiment. 
Following Ref. ||^ we choose (3 = 2. In Ref. this choice was motivated by the phase of the 
quantum oscillations. Using the approach of Fleischmann et al. ||2^ we have also computed 



numerically the classical resistivity for this potential. Experimentally one finds a single peak 
located at Rc = a/2. As shown in Fig. 4, the model potential with (3 = 2 also leads to a 
single dominant peak at the same magnetic-field value in the classical resistivity. This lends 
additional support to this choice of (3. Other choices such as (3 = A lead to considerable 
additional structure in the resistivity |^ . 



To study the nature of the classical dynamics for this choice of model potential we have 
generated Poincare surfaces of section. For (3 = 2 and magnetic fields such that Rc = a/2 
one finds that phase space is entirely chaotic (within our resolution). This is shown in 
Fig. 5(a). The classical dynamics becomes more (less) chaotic as one decreases (increases) 
the magnetic field from Rc = a/2. In particular, islands of regular motion appear for 
larger magnetic fields. Our results for completely chaotic systems should be applicable for 
magnetic fields with Rc > a/2. However, as already mentioned above, one might expect 
a somewhat reduced level of accuracy because the trajectories are close to stability. The 
classical dynamics also becomes less chaotic as one increases the steepness of the antidot 
potential. This is shown in Fig. 5(b) for (3 = 4 and Rc = a/2, where one finds a large regular 
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island in the Poincare section. 

Our central result Eq. (|36| ) expresses the quantum oscillations in the conductivity in 
terms of the classical periodic orbits of the antidot lattice. Since long orbits are suppressed 
due to elastic and inelastic scattering, stability amplitude, and finite temperature, we focus 
on the contributions from the shortest periodic orbits labeled by (a), (b), and (c) in Fig. 6. 
It turns out that the action of orbit (a) changes only slowly with magnetic field giving rise 
to oscillations with large period. These are not resolved experimentally because of the small 
magnetic-field range. Note that while orbit (b) has fourfold symmetry, orbit (c) has only 
twofold symmetry so that the latter contributes twice. As already pointed out by Weiss et 
al. these orbits explain the phase of the quantum oscillations. Both orbits have nearly 
the same action and one finds that the predicted minima from Eq. (|36|) indicated by arrows 
in Fig. 3 agree well with the observed minima in the conductivity. Furthermore, our results 
predict in agreement with experiment that longitudinal and Hall conductivity are in phase. 

The amplitude of each periodic orbit 7 in Eq. ( p6D depends on the traversal time Ty, the 
velocity correlations C^{vx, v^), and the stability factor |TrMy — 2|^^/^. Our numerical results 
for these quantities at the classical peak = a/2 and for uJcT^i = 2 are collected in table 
1. Orbit (c) in fact consists of three periodic orbits in close vicinity of each other. In view 
of the Fermi wavelength Ai? we do not believe that it is justified to include all three orbits 
separately. Instead we only include one of the three with typical stability. Presumably 
this tends to underestimate the amplitude and hence a more accurate treatment would 
lead to closer agreement between theory and experiment: Using the numbers of table I 
we find 6a^^ ^ 0.008(kfi)-i from Eqs. (|6D and 6(T^y ^ 0.004(k(])-^ from Eq. (||) |8 



These theoretical estimates are smaller by factors of two and four, respectively, than the 
experimental results. We believe that this order-of- magnitude agreement is reasonable given 
the theoretical uncertainties mentioned at the beginning of this section. 

Our results predict a beat period due to the Zeeman splitting which is rather large 
compared to the magnetic-field range over which the quantum oscillations were observed. 
While the effects of the Zeeman splitting were thus unobservable, the predicted temperature 
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dependence is consistent with experimental results. The amplitude of the quantum oscil- 
lations is expected to decrease exponentially with temperature on the scale K/tiT^. Hence, 
one predicts a characteristic temperature of 0.7K. This is consistent with experiment where 
quantum oscillations were observed at 0.4K but not at 4.7K [^Q. 

In the introduction we remarked that the periodicity of the quantum oscillations varies 
between 1/B periodicity and B periodicity depending on the experimental sample DH- To 
understand this observation, consider the magnetic-field range where the cyclotron radius is 
of the order of half the lattice spacing so that electrons can circle around a single antidot. 
If the antidot diameter is small the relevant periodic orbits are only weakly perturbed by 
the presence of the antidot potential. Hence, one expects 1/B periodicity familiar from the 
Shubnikov-de Haas oscillations in unpatterned samples. On the other hand, if the antidot 
diameter is larger, the periodic orbits are modified by the antidot potential. In particular, 
they will be inhibited in contracting with increasing magnetic field. Hence, with increasing 
antidot diameter one crosses over to a situation where the enclosed area of the periodic 
orbits is approximately constant with magnetic field. In this limit the change of the classical 
action with magnetic field is merely due to the Aharonov-Bohm flux and hence one finds B 
periodicity. 

IV. DISCUSSION AND SUMMARY 

Motivated by a recent experiment by Weiss et al. ||^ we have studied the quantum 
transport in large antidot lattices. We have evaluated the Kubo formula for the longitudinal 
and the Hall conductivity semi classically for systems whose classical dynamics is chaotic. 
While the leading contribution of the semiclassical expansion is just the classical Kubo 
formula as expected, we have identified quantum corrections due to the classical periodic 
orbits of the system. The contribution from each periodic orbit oscillates as function of 
the Fermi energy and the magnetic field with the phase given by the classical action of the 
orbit. The amplitude is determined by the stability and the velocity correlations of the orbit. 
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Furthermore, we find that the periodic-orbit contributions to the conductivity decrease with 
increasing temperature on the scale h/T^ where denotes the traversal time of the orbit. 
Our results provide a consistent explanation for the quantum oscillations superimposed on 
the classical magnetoresistivity of antidot lattices recently observed by Weiss et al. 

Our approach implies that the quantum oscillations in antidot lattices are closely related 
to the well-known Shubnikov-de Haas oscillations. We have made this relation explicit 
by giving an analogous semiclassical treatment of Shubnikov-de Haas oscillations. One 
finds that the amplitude of Shubnikov-de Haas oscillations is of a lower power in h and 
hence larger than that of the quantum oscillations in chaotic systems. In this sense, the 
periodic-orbit contribution to the conductivity is a sensitive probe of quantum chaos. This 
is analogous to the well-known result that periodic-orbit contributions to the density of states 
are much larger for integrable systems (Berry- Tabor formula |^) than for chaotic systems 
(Gutzwiller formula |2^). Most importantly, the close relationship between Shubnikov-de 
Haas oscillations and the periodic-orbit contribution in chaotic systems implies that both 
effects should be observable under analogous conditions. Hence observation of the latter 
requires that the elastic and inelastic mean free paths be much larger than the length of the 
periodic orbit but may be smaller than the sample size. Furthermore, the temperature must 
be smaller than h/T^ which is the analog of the Landau-level spacing for Shubnikov-de Haas 
oscillations. Both effects also require a macroscopically homogeneous sample. Otherwise 
periodic orbits in different parts of the sample contribute with different phases resulting in 
a reduced amplitude due to destructive interference. 

Our calculation of the periodic-orbit contribution to the conductivity of antidot lattices 
follows the spirit of but is different in detail from Gutzwiller's derivation of the periodic- 
orbit contributions to the density of states. Like Gutzwiller's trace formula, our results are 
strictly valid only for completely chaotic systems. Of course, many systems exhibit mixed 
dynamics with both regular and chaotic regions in phase space. In fact, antidot lattices 
tend to become mixed systems with increasing magnetic field or increasing steepness of 
the antidot potential. A quantitative theory for the periodic-orbit contributions to the 

24 



conductivity in mixed systems remains a challenge for future work. Fortunately, the antidot 
sample which was used to demonstrate the quantum oscillations for the first time [^] appears 
to be quite close to completely chaotic dynamics over a significant magnetic-field range so 
that a comparison with our theoretical results is meaningful. 

It is interesting to compare our work with other recent semiclassical studies of transport 
properties. Physically, one also expects a weak-localization correction to the classical con- 
ductivity for the antidot lattice. It appears reasonable that the weak-localization correction 
for antidot lattices should be of the same order in Planck's constant h as that for diffusive 
systems. This would imply that the weak-localization correction is of the same order in h as 
the periodic-orbit contributions computed in this paper. Nevertheless, the stationary-phase 
conditions do not allow for interference between time-reversed trajectories. Recently, Arga- 
man argued in the context of diffusive systems that a resolution of this dilemma requires 
one to go beyond standard semiclassical theory. For small but finite h he also included 
trajectories which approximately (but not strictly) satisfy the stationary-phase conditions 
and was then able to rederive the diagrammatic results for weak localization in the diffusive 
regime. 

Semiclassical expressions for the conductance have also been derived recently starting 
from the Landauer-Biittiker formula for the two-probe conductance. In particular, semi- 
classical treatments have been given for weak localization and conductance fluctuations in 
phase-coherent ballistic microstructures. In this formulation the semiclassical conductance 
is given in terms of a double sum over classical trajectories entering through one lead and 
exiting through the other. It is not clear how periodic orbits of the microstructure would 
enter in this approach. It is clearly desirable to gain a deeper understanding of the relation- 
ship between these different approaches to judge conclusively the range of validity of the 
semiclassical results for transport properties. 
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APPENDIX A: SEMICLASSICAL GREEN FUNCTION FOR ELECTRONS IN A 

MAGNETIC FIELD 

In this appendix, we derive the semiclassical Green function in Eq. (^) for an electron 
in the two-dimensional x-y plane subject to a perpendicular magnetic field. We include a 
weak disorder potential resulting in the broadening F = h/2Tei of the Green function defined 
above Eq. (^. We recall that the Green function is the Fourier transform of the propagator, 

Gt'AE) = -J ^dtexp + tr)tj K{r', r; t), (Al) 

where the propagator 

K{r',r;t) = e{t) {r'\exp{-iHt/n)\r), (A2) 

is known exactly j41||. 



r — r'P Urt 



cot {xy' — x'y) 



2 2 



(A3) 



To simplify notation we denote the phase of if(r', r; t) by 

R{r, r'; t) = - f ' ^ ' cot -j- - {xy' - x'y)\ . (A4) 

In the semiclassical limit ^ ^ the integral in Eq. ( [Al| ) can be performed by the stationary- 
phase method. The stationary-phase condition reads 

j^_mujl( |r-r'| ^ 

~ ~8~ Vsin(tUct/2) ) ■ ^ ' 
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This condition has no solutions for |r — r'| > 2i?c reflecting the fact that there are no 
classical trajectories between two points further apart than twice the cyclotron radius Re- 
Thus the semiclassical G^, t.{E) vanishes for |r — r'| > 2Rc. On the other hand, points with 
|r — r'l < 2Rc lie on two different cyclotron orbits (cf. Fig. 1), and one finds an infinite set 
of stationary times Tn^q, n = 0,1,2, q = S, L, 



T 



2Tm 



n,q 



1 ' 



(A6) 



with 



ts = — arcsin(|r — y'\/2Rc), 



t, 



TT — arcsm r 



UJr 



r'|/2i?e)]. 



(AT) 



Here n counts the number of complete cyclotron revolutions and q distinguishes between the 
two orbits. The T^^q are the traversal times for the classical motion from r to r'. Performing 
the stationary-phase integrals, one obtains 

1 \ V- - ^ / d^R 



Gt,JE) 



m 



2ih \27iih 



1/2 oo 

h^qhL sin(^cr„,,/2) \ 



-1/2 



X exp(-r„^g/2re;) exp |^ \ETr,^q + i?(r, r'; T„„g)]| . 



One has 



d'^R 



uJrE cot 



(A8) 



(A9) 



and replaces all trigonometric functions in the amplitude by their absolute values. The 
phases are absorbed into the Maslov indices 



The phase appearing in Eq. 



5(r,r';^) =ET„,,+ 



rin,s = 2n, 
r]n,L = 2n + l. 
is just the action 



■ cot 



- {xy' - x'y) 



where the stationary times T„^g depend on r, r', and E as described by Eqs. 
everything together one arrives at the result given in Eq. (P). 



(AlO) 



(All) 



Putting 
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APPENDIX B: STABILITY AMPLITUDE 



In this appendix a derivation of Eq. (p^ is presented. First, we note that evaluation 
of the amphtudes of the Green function in the local coordinate system [y, z) of the 
trajectory yields [El] 



dydy' 



(Bl) 



Hence we need to simplify the expression 



71 



dydy' 



72 



dydy' 



1/2 



detFl 



(B2) 



To do so, we introduce a matrix A^*) for trajectory 7j which relates to linear order small 
transverse deviations in the initial and final momentum = [6p, 6p'] to small transverse 
deviations in the initial and final positions Sy'^ = [6y, 6y'] , 



5p = A5y. 



(B3) 



In terms of the action S^^ the matrix A*^*^ can be readily expressed as 

d^S^ 



d^S. 



_ d'S,^ 
dydy' 
d^S^^ 



(B4) 



dy'dy dy''^ 

Hence, all quantities appearing in a-^ can be expressed in terms of matrix elements of A^^^ 
and A^'^\ It turns out that is more conveniently expressed in terms of the monodromy 
matrices M and of the trajectories 71 and 72, respectively. The monodromy matrix 
relates to linear order the final transverse deviations (5x'^ = [5y' , 5p'] to the initial transverse 
deviations in phase space 5x-^ = [5y, 5p\ , 

5x' = M5x. (B5) 



An important property of the monodromy matrix is that it has unit determinant, det M = '. 
The matrices A^^^ and A^'^'' can now be easily written in terms of the matrix elements m. 
and Uij of the monodromy matrices, respectively. One finds 



28 



(B6) 

and likewise for A^'^\ Using these relations we find for in terms of the monodromy 

matrices M and A^, 

(l/mi2)(lM2) 

(mii/mi2 - nii/ni2){m22/mi2 - ^22/^12) - (l/"^i2 - 1/^12)^ 

nl2{mum22 - 1) + mi2{nun22 - 1) - ni2mi2{mun22 + ^11^22 - 2) 
'n'i2"^2i + '^i2'n'2i - "^11^22 - niim22 + 2|"^ 

TtMN'^ - 2\-\ (B7) 

Here we used that det M = det = 1. Finally, we note that MN~'^ = M^, where is the 
monodromy matrix of the (primitive) periodic orbit. Hence, one finds the result 

a-y = \TtMP - 2\-^/\ (B8) 

which proves Eq. (p4D. It is remarkable that only the stability of the periodic orbit enters 
in this factor, regardless of how long the trajectories 71 and 72 are. 



^(1) 



mil 
mi2 



mi2 



m2i- 



miim22 m22 
mi2 mi2 




APPENDIX C: TERMS INVOLVING TWO ADVANCED (RETARDED) GREEN 

FUNCTIONS 

In this appendix it is shown that the contributions to the conductivity (0) involving two 
advanced or two retarded Green functions vanish in the semiclassical approximation. While 
we focus on the calculation for antidot lattices, the same statement holds in our treatment of 
Shubnikov-de Haas oscillations because the factor displayed in ( |C2|) below is also present in 



this case. For definiteness, consider the contribution a^J^. Using the semiclassical expression 
for the Green function (|19[) one has 



(r',r) 

exp [jrlS^^i'^, r') + S^,{r', r)] - y [r/^^ + ^71]} • (CI) 



7l(r,r') 72(r',r) 
X 
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One notes that in contrast to this contribution involves both trajectories running from 
r to r' and from r' to r. In particular, there is no classical contribution for which the phase 
factor cancels. Again performing the spatial integrals by the stationary-phase method, one 
finds that the integrand becomes stationary when vi = V2 and = Vg, where Vj and 
denote the velocities of 7^ at position r and r', respectively. The stationary-phase condition is 
satisfied when both 71 and 72 are parts of the same periodic orbit. This is illustrated in Fig. 
7. A straight-forward extension of the calculation in appendix shows that the amplitude 
factors also combine in this case to yield |TrM^ — 2|^^/^ involving only the monodromy 
matrix of the periodic orbit. However, the contribution a^^ vanishes in the semiclassical 
limit because of the integrals over z and z' along the periodic orbit. The analog of Eq. (|35D 
reads 

00 .'IP .'IP 

exp(-pr^/2rez) j j dt' v^{t)v^{t + t') exp{-nTjTei) (C2) 

This expression vanishes because the exponential suppression factor due to disorder does 
not depend on t or t' and hence the integrals are over one period of a periodic function with 
zero mean. A completely analogous argument can be used to show that a~~ vanishes in the 
semiclassical approximation. 
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FIGURES 



Fig. 1 Sketch of classical trajectories which contribute to the semiclassical Green function 
G^, r(-E') for electrons in a homogeneous magnetic field. The two shortest trajectories 
from r to r' labeled by S (for short) and L (for long) are depicted by bold lines. There 
are also contributions from trajectories which complete additional revolutions around 
the cyclotron orbits of radius Re- 

Fig. 2 Sketch of classical trajectories 71 and 72 from r to r' which satisfy the stationary- 



phase conditions (|29| , pOl) . Both trajectories follow a periodic orbit 7 (dashed line). 



The shaded depicts an antidot. 

Fig. 3 Experimental results of Weiss et al. |^,37] for quantum oscillations a{T = 0.4K) — 
cr(T = 4.7K) as function of magnetic field B in longitudinal conductivity axx (full line) 
and Hall conductivity axy (dashed line). The arrows indicate the magnetic- field values 
where minima in the conductivity are predicted according to Eq. (BBf). 



Fig. 4 Numerical result for classical resistivity of the model potential ( p^ with jS = 2 
as function of magnetic field B. The resistivity exhibits a single dominant peak at 
B = Bq where the cyclotron orbit equals half the lattice spacing of the antidot arrays. 
The parameters were chosen according to the experimental parameters for the main 
sample investigated in Ref. which also showed a single peak in the resistivity at the 
same magnetic-field value. This supports the choice of /5 = 2 for the steepness of the 
antidot potential. 



Fig. 5 (a) Poincare surface of section [x/ a,Vx/ \j Ep/m) at ?/ = (mod a) with exponent 
(3 = 2 for i?c = 0/2 and antidot diameter d = a/2 appropriate for the experimental 
sample in Ref. 0. Within our resolution the classical dynamics is completely chaotic, 
(b) The same surface of section for (3 = 4: and d = aj?). Here one finds a mixed phase 
space with both islands of regular motion and regions of chaotic motion. 
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Fig. 6 Contour plot of the antidot potential (thin hnes) together with the relevant periodic 
orbits for B = Bq where the cyclotron radius Rc = a/2. The innermost lines correspond 
to the boundaries of the classically-forbidden regions. 

Fig. 7 Sketch of classical trajectories 71 from r' to r and 72 from r to r' which satisfy 
the stationary-phase conditions for the terms involving two advanced or two retarded 
Green functions. Both trajectories follow a periodic orbit 7 (dashed line). The shaded 
area depicts an antidot. 
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TABLES 



Table I Numerical results for the traversal times T^, velocity correlations Cy{vx,Vx), and 
stabilities |TrMy — 2|^^/^ of the periodic orbits shown in Fig. 6. The parameters 
are chosen appropriately for the experimental sample in Ref. The magnetic field 
corresponds to Rc = a/2. For comparison, numbers are also given for the cyclotron 
orbits of unpatterned samples. 



Orbit 






|TrM^ - 2|-i/2 


b 


3.42 


1.41 


0.42 


Cl 


4.04 


1.62 


0.62 


C2 


4.07 


1.63 


0.46 


C3 


4.04 


1.62 


0.62 


cyclotron 


3,14 


1.26 




orbit 
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